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ABSTRACT 

We have worked out predictions for the radio counts of star-forming galaxies down to nJy levels, 
along with redshift distributions down to the detection limits of the phase 1 Square Kilometer Array 
MID telescope (SKA1-MID) and of its precursors. Such predictions were obtained by coupling epoch 
dependent star formation rate (SFR) functions with relations between SFR and radio (synchrotron and 
free-free) emission. The SFR functions were derived taking into account both the dust obscured and 
the unobscured star-formation, by combining far-infrared (FIR), ultra-violet (UV) and Ha luminosity 
functions up to high redshifts. We have also revisited the South Pole Telescope (SPT) counts of dusty 
galaxies at 95 GHz performing a detailed analysis of the Spectral Energy Distributions (SEDs). Our 
results show that the deepest SKA1-MID surveys will detect high -2 galaxies with SFRs two orders 
of magnitude lower compared to Herschel surveys. The highest redshift tails of the distributions at 
the detection limits of planned SKA1-MID surveys comprise a substantial fraction of strongly lensed 
galaxies. We predict that a survey down to 0.25 pJy at 1.4 GHz will detect about 1200 strongly 
lensed galaxies per square degree, at redshifts of up to 10. For about 30% of them the SKA1-MID 
will detect at least 2 images. The SKA1-MID will thus provide a comprehensive view of the star 
formation history throughout the re-ionization epoch, unaffected by dust extinction. We have also 
provided specific predictions for the EMU/ASKAP and MIGHTEE/MeerKAT surveys. 

Keywords: radiation mechanisms: general - radio continuum: galaxies - galaxies: evolution - galaxies: 
starburst - infrared: galaxies 


1. INTRODUCTION 

There are several well established diagn ostic tools for measu ring the star formation rates (SFRs) in galaxies, each 
with its own strengths and weaknesses (see iKennicutt fc Evausll2012L for a comprehensive review). The UV and the Ha 
luminosities are among the best direct tracers of the emission from young stars, so that the SFR is simply proportional 
to the luminosity. These diagnostics, however, are highly sensitive to dust extinction and may miss obscured star 
formation. The latter is measured by the far-infrared/sub-mm luminosity due to dust reprocessing of the emission 
from newly formed stars. But the dust heating is not always dominated by young stars: older stellar populations can 
contribute as well. Both the UV/Ha and the far-infrared (FIR) diagnostics may also be contaminated by emissions 
from Active Galactic Nuclei (AGNs). _ 

Ever since radio surveys have reached sub-mJy flux density levels dWindhor st et al.1 ll985T) . they proved to be a 
primary means of identifying star-forming galaxies at high redshift ('e.g.. IBonzini et alj|2013f) . Observations have 
solidly demonstrated a tight correlation between the low-frequency radio continuum, due to synchrotron emission of 
relativistic electro ns mostly produced by supernovae, and f ar-IR emission of galaxies, legitimating its application as 
a SFR tracer (see iHelou. Soifer. fe Rowan- Robinsonl 119851: iCondonl [19921 : lYun. Reddy fc Condonl 200ll : llvison et al.l 
120101 IJarvis et al.l I2010L IBourne et al.l 1201 ll IMao et al.l 1201 If) . A word of caution is in order, however, because the 
physical basis of this relation is not fully understood. Actually there are many physical processes, such as propagation 
of relativistic electrons, strength a nd st ructure of the mag netic field, size and composition of dust grains, that must 
cooperate to produce this relation (lBelIll2003l : lLacki et al.ll2QToh . As a consequence it is not granted that this relation 
applies also to redshift/luminosity ranges where the available data are insufficient to test it precisely. In addition the 
synchrotron emission could be contaminated, even by a large factor, by faint nuclear radio activity, and this would 
bias the estimation of the SFR. 

A more direct radio SFR tracer is the free-free emission from the gas ionized by massive young stars (lMurphvll2009i ; 
iMurphv et al.il2014 b Being directly proportional to the production rate of ionising photons it provides a measure of 
the SFR without the complication of dust attenuation encountered in the optical/UV. It has a flat spectrum and is 
expected to show up at frequencies of tens of G Hz. The interpretatio n of d ata at these frequencies may be complicated 
by the presence of “ano malous” dust emission dPlanck Collaboration XXH201H . and references therein) attributed to 
spinning dust grains (e.g. JDraine fc Lazarianlll998f b However a significant contribution of this component to the global 
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emission of galaxies has not been proved yet dMur phv et aDl2012t iPlanck Collaboration XXVll2014D . 

The current deepest radio surveys Csee lde Zotti et al.ll2010t for a review) have not been carried out at frequencies 
high enough or are not deep enough to see the transition from the synchrotron to the free-free dominance. Only with 
the advent of the Square Kilometer Array (SKA), and partially with its precursors, we expect to efficiently select 
high -2 star-forming galaxies via their free-free emission. 

In this paper we carry out a thorough investigation of the radio counts and redshift distributions of star forming 
gala xies, down to the flux density levels that can be reached by the SKA. Our predictions rely on a state-of-the-art model 
(iCai et al.ll2014h that accounts for the evolution of the cosmic SFR function up to z ~ 10 as measured by a combination 
of far-infrared (FIR), UV and Ha data. An outline of the model is presented in Sect. [2] The redshift dependent SFR 
function is translated into the evolving radio luminosity function, taking into account both synchrotron and free-free 
emissions, using the calibrations bv [Murphv et al. (201 If ). In Sect. Owe tune the correlation between the SFR and 
the synchrotron lum inosity in order to match the observational determination of the local luminosity function of star¬ 
forming galaxies (|Mauch fe Sadled 1200 71 . This allows us to reproduce the available multi-fre quency source counts . 
Those of star-forming galaxies at 95 GHz provided by the South Pole Telescope (SPT) survey (|Mocanu et al.l feOlS). 
that are the most sensitive to the free-free emission, have been reassessed performing a more accurate discrimination 
between star-forming galaxies and radio-loud active galactic nuclei. Having the radio luminosity function at each 
redshift it is straightforward to work out predictions for counts and redshift distributions of star-forming galaxies. In 
Sect. [4] we show a selection of the results for the frequency range 1.4 - 30 GHz and compare the coverage of the SFR -2 
plane by Herschel, optical/UV surveys and the phase 1 SKA-MID. Section [5] explores possible synergies with deep 
optical/near-infrared (near-IR) surveys. Finally, Sect. [ 6 ] summarizes our main conclusions. 

Throughout this paper we adopt a flat ACDM cosmology with matter density fi m = 0.32, fib = 0.049, Ha = 0.68, 
Hubble constant h = Ho/lOOki ns -1 M pc -1 = 0.67 , spec trum of primordial perturbations with index n = 0.96 and 
normalization erg = 0.83 (Planck Gollaboration XVli20l4 ). 

2. OUTLINE OF THE MODEL 

As illustrated by the right-hand panel of Fig. 8 of iMadau fc Dickinsonl (2014f) and by Fig. 10 of iGai et~ahl (20141) . 
the mean dust attenuation is low at z > 6 , when the metallicity of the intergalactic medium (IGM) is lower than the 
solar value, Z@, by a factor < 10 -3 (Fig. 14 of IMadau fo Dickinsonl [201 41. Thus, at these redsliifts, measurements at 
UV wavelengths are likely to be a reliable SFR di agnostic for most ga l axies, although a fe w galaxies with high SFRs 
and substantial dust extinction have been found (Riechcrs et al.1 120131 : 1 Watson et al.ll2015D . Gr eat progress has been 
made in recent years in the determination of UV luminosity functions at high 2 , up to 2 ~ 10 (Bouwens et aI1l2008L 
l2011afbL 120151 : iSmit et aDl2012t lOesch et al.ll 2012 t I2013L12014 iSchenker et al.ll2013l : iMcLure et al.ll2013D . struggling to 
reconstruct the history of the cosmic re-ionization. 

The mean dust attenuation increases rapidly with decreasing 2 down to z ~ 1—1.5, where it peaks, to decline at 
still lower 2 . In the range 1 < 2 < 3, including the peak in the cosmic SFR density, the star formation is heavily 
dust enshrouded, so that the IR ( 8 - 1000 /im) luminosity is an accurate SFR measure. At the other redshifts the safest 
approach combines UV with IR measurements. In dealing with the latter, we need to be aware that, especially at low 
redshifts, a significant contribution to dust heating may be due to older stellar populations. 

A detailed study of the evolution of the SFR function across the cosmic history was carried out bv lCai et all (|2013f ) 
focussing on IR data, and by ICai et all (|2014D at 2 > 2 focussing on U V and Lyq data but t aking into account also 
dust atten u ation and re-emis si on. These studie s build on the work bv iGranato et al.1 ()2004D . further developed by 
lLapi et, all (2006112011L120141) . iMao et all (20071). The unde r lying scenario relies on the results of inte nsive IV-body 
simulations and semi-analytic studies (e.g.. iZhao et al.l l2003t iWang et a.l.ll 2011 t lLapi fc Cava.lierdl 20 TTI) showing that 
pre-galactic halos undergo an early phase of fast collapse, including a few major merger events, during which the 
central regions reach rapidly a dynamical quasi-equilibrium. A slow growth of the halo outskirts in the form of many 
minor mergers and diffuse accretion follows. This second stage has little effect on the inner potential well where the 
visible galaxy resides. 

The baryons falling into the potential wells created during the fast collapse phase are shock-heated to the virial 
temperature. The gas cools down giving rise to star formation. At the same time it loses angular momentum, e.g. 
by effect of the radiation drag, and flows towards the central regions, accreting onto the centr a l super-ma ssive black 
hole. These processes are governed by a set of equations, reported in the Appendix of ICai et al.1 (2013D . that take 
into account the energy feedback from supernovae and from the active nucleus. The model includes a self consistent 
treatment of the chemical evolution of the IS M, cal cul ated us ing the standard equations and stellar nucleosynthesis 
prescriptions, as written down, for instance, in iRomano et all (2.002). The chemical evolution controls the evolution 
of the dust abunda nce, hence the dust absorption and re-emission. 

In this framework ICai et al.l (|2014t ) worked out a model for the evolution of the ultraviolet (UV) and Lya luminosity 
functions of high-redshift (2 > 2 ) galaxies, taking into account their chemical evolution and the associated evolution 
of dust extinction. The model yields good fits of the observed luminosity functions at all redshifts at which they 
have been measured and provides a physical explanation for the weak observed evolution of both luminosity functions 
between 2 = 2 and 2 = 6 . 

While at high 2 the high star formation efficiency impl ied by the s cenario briefly sketched above applies to most 
galaxies, at lower 2 a dichotomy arises (see, e.g.. iBernardi et aLll 2010 D . Early-type galaxies and massive bulges of Sa 
galaxies are composed of relatively old stellar populations with mass-weighted ages > 8-9 Gyr. They must therefore 
have been forming most of their stars over short timescales at 2 > 1 1.5 and must have dominated the cosmic SFR 
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at high z. Instead, the disc components of spirals and the irregular galaxies are characterized by significantly younger 
stellar populations and star formation activity continuing up to the present time. They are therefore the dominant 
star formers at z < 1.5. 

iCai et all 1 2013 1 modeled the cosmological evolution of the SFR function of proto-spheroidal galaxies (the progenitors 
of present day early-type galaxies and of massive bulges of disc galaxies) during their dust-enshrouded star-formation 
phase, and the concomitant evolution of their active nuclei. This model provides a physical explanation for the 
observed positive evolution of both galaxies and AGNs up to z ~ 2.5 and for the negative evolution at higher redshifts, 
for the sharp transition from Euclidean to extremely steep counts at (sub-)mm wavelengths, as well as for the (sub- 


approaches. Furthermore, as shown bvIXia et ah (j2012f) and ICai et al. 

(2013 

). the halo masses inferred both fro 

m the 
power 

angular correlation function of detected sub-mm galaxies (Coorav et al 

20101 

Maddox et al.||2010D and from the 

spectrum of fluctuations of the cosmic infrared background (Amblard et al.l 

2011: Planck Collaboration XVIII 

20 lH 


are larger than th ose implied by the major mergers plus top-heavy in itial stellar mass function (Baugh et alJl200a 
lLacev et al. 201$ and smaller than those implied by cold flow models dDave et alJl2010h . 

Late-type galax ies fo llow a differ ent evolutionary track. The evolution of the IR luminosity function of these objects 
was described by ICai et ahl ( 2013 ) using a phenomenological approach, considering two populations with different 
spectral energy distributions (SEDs) and different evolutionary properties: “normal” late-type galaxies, with low 
evolution and low dust temperatures and rapidly evolving starburst galaxies, with warmer dust temperatures. 

The evolution of late-type “normal” and “starburst” galaxies is described using a phenomenological parametric 
approach. For the starburst population, power-law density and luminosity evolution [$*( 2 ) = <l>*(z = 0) x (1 + z ) a *; 
A*(z) = L*(z = 0) x (l+z) aL ] up to ^break = 1 was assumed. The normal galaxy population has only a mild luminosity 
evolution up to the same value of Zbreak, as indicated by chemo/spectrophotometric evolution models. At z > Zbreak 
both <f>*(z) and A*(z) are kept to t he value s at Zh rsak multiplied by a smooth cut-off function. For further details and 
the values of the parameters, see ICai et al.1 (120131) . 

The model fits a broad variety of data 5 : multi-frequency and multi-epoch luminosity functions of galaxies and AGNs, 
redshift distributions, number counts (total and per redshift bins). It also accurately accounts for t he recent counts 
and redshift dist ribution of st rongly lensed galaxie s detected by the South Pole Telescope 1SPT: IMocanu et ahll2013l : 
I Weill et aLll2013T) as shown bv iBonato et al.l (|2014D . 

When the dust heating is dominated by young stars and the effective optical/UV optical depth is high, the total 
(8-1000/tm) IR luminosity, Air, is a good measure of the SFR. More generally, the IR luminosity is the sum of a 
“warm” dust component heated by young stars and of a “cold” dust (or “cirrus”) component, heated by the general 
radiation field probably domin ated b y olde r st ars. Obviously the only tracer of the SFR is the first component. This 
issue was investigated bv lClemens et al.l (120131) using a complete sample of local star-forming galaxies detected by the 
Planck satellite. We adopt the relation between SFR and Air derived by these authors: 


log(SFR/M 0 yr *) = -9.6 + log(L IR /L 0 ) 




2.0 


Vlog(Li R /L Q ) - 7.0 


(1) 


As expected, high IR luminosities generally correspond to high SFRs, i.e. to a dominant “warm” dust component. Note 
that calibrations of SFR dia gnostics g e nerall y depend on the adopte d initial mass function (IMF ). In the above equation 
we have adopted the one byl Chabricr (2003 ). According to Fig. 1 oflCho miuk fe Povichl (1 20111) the differences with the 
calibrations usi ng thelKroutia J2001I) IMF ( Murphy et al.ll201ll 120121 : iKennicutt fc Evansll2012lf are < 10% and will be 
neglected. The ICai et al.l (I2013H model deals only with dust-obscured star formation. To get a comprehensive view of 
the evolution of the SFR function of late-type galaxies, including the unobscured star formation, we have complemented 
it taking into account optical/UV SFR indicators, adopting again a parametric approach for evolution 6 . Details are 
given in the Appendix [A] 

The SFR functions at several redshifts, obtained combining IR, UV and Ha su rvey data, are shown in Fig. [TJ where 
also the SFR functions derived from the UV luminosity functions using the IKennicutt fe Evans! (l2012h calibration but 
without the extinction correction are plotted (thin solid lines). Below a SFR ~ 0.1-1 Mq yr~ , depending on z, the 
SFR functions derived from the IR luminosity functions fall below those derived from the UV and Ha luminosity 
functions without extinction corrections, implying that extinction becomes negligible for low SFRs. The opposite 
happens for higher SFRs. The model tends to underpredict the brightest points of the luminosity functions. We note, 
however, that the discrepancy is of the order of differences among estimates from different surveys. This suggests that 
the data points are affected by cosmic variance due to the relatively small area studied; the rare objects with the 
highest luminosities are preferentially detected in areas where they are overdense. 

While at z < 1.5 the observationally determined SFR functions are fully accounted for by UV measurements with 
the appropriate extinction corrections, at higher redshifts UV surveys miss most of the galaxies with the highest SFRs, 
that are heavily dust obscured and therefore show up at far-IR/sub-mm wavelengths. According to our model, the 
dominant contribution at the highest SFRs at z > 1.5 comes from proto-spheroidal galaxies. 

5 See http://people.sissa.it/~zcai/galaxy_agn/ or http://staff.ustc.edu.cn/~zcai/galaxy_agn/index.html 

6 Early-type galaxies are in passive evolution at low z and their early unobscured star-formation phases are dealt with by ICai et ahl 
(l20T4h 
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Figure [2] presents a synoptic view of the cosmological evolution of the cosmic SFR function up to z = 8 as given by 
the model described above, including the contributions of normal late-type, starburst and proto-spheroidal galaxies. 
This function has a complex evolutionary behaviour. Its knee value increases w_ith_z up to z ~ 2- 3 and declines at 
higher redshifts; its slopes above and below the knee also change with redshift. lAversa et al.1 ( 2015 1 have worked out 
approximate analytic expressions providing a sensible rendition of the observational determinations of the redshift- 
dependent SFR function. To this end they used a modified Scliechter function with evolving characteristic luminosity 
and slopes. The luminosity function [Mpc -3 log(SFR) -1 ] writes 


where 


</>SFR (z) = $(z) 


SFR 
SFR Jz) 


-| l-a(z) 


exp 


SFR 


Wz)' 


SFR*(z)_ 


log<£( 2 ) [Mpc 3 log(SFR) x ] = —2.4 — 2.3% + 6.2 x 2 — 4.9% 3 , 
log SFR* ( 2 ) [Mr., yi- 1 ] = 1.1 + 3.2 * - 1.4 X 2 - 2.1 X 3 , 
a ( z ) = 1-2 + 0.5 x — 0.5 x 2 + 0.2 x 3 , 
cj(z) = 0.7 - 0.15 x + 0.16 x 2 + 0.01 y 3 , 


( 2 ) 


with x = log(l + z). Note that the above equation is simply meant to provide an analy tic renditio n of th e complex 
evolution of the SFR function, for use in practical applications. As pointed out by lAversa et al.l (|2015f l the non- 
homogeneous nature of the data sets used to derive the SFR functions and the diverse systematics affecting them do 
not warrant the use of formal fitting procedures. 


3. RELATIONS BETWEEN RADIO EMISSION AND SFR 
3.1. Calibration of the relation between SFR and synchrotron emission 

A calibration of the SFR-s ynchrotron lum inosity relation was obtained by [Murphy et al.l (|2011l) using Starburst99 
IjLeitherer et alJll999f) with a lKrouoal (1200111 IMF and specific choices for the metallicity and the star formation history. 
We have slightly modified their relation including a steepening of t he synch rotro n spect rum by Aa = 0.5 above a break 
frequency of 20 GHz to take into account electron ageing effects (jBandav fe WolfendalellT99ll) . The SFR-synchrotron 
luminosity relation then writes: 


L sync ~1.9x 10 28 


/ SFR \ 

Uioyi-V 



erg s 1 Hz 1 


(3) 


The high-frequency synchrotron emission can be further suppressed at high 2 because of the energy losses of relativistic 
electrons due to inverse Compton scattering off the Cosmic Microwave Background (CMB) photons, whose energy 
density increases as (1 + 2) 4 (jNorris et al.ll2013L iCarilli et al.ll2008t lMurphvll2009H . Also it should be noted that, at 
variance with free-free and UV emissions which are a measure of the current SFR, the synchrotron emission is delayed 
by ~ 10 7 yr, the main sequence lifetime of massive stars whose supernova explosions are responsible for the production 
of the relativistic electrons (|Clemens et alJ|2008D . At high 2 the fraction of galaxies younger than ~ 10' yr, which are 
lacking relativistic electrons, hence the synchrotron emission, increases. 

Since these effects are not taken into account in our calculations because we lack an adequate model for them, we 
may overestimate the faint counts and the high -2 tails of the redshift distributions. The overestimate of the total radio 
luminosity is however unlikely to be large because even at 1.4 GHz (and more so at higher frequencies) at flux densities 
below the current observational limits the free-free contribution is comparable to the synchrotron one (see Fig.[3|. 

In translating the SFR function into a synchrotron luminosity function we need to take into account that the corre¬ 
lation between the SFR and the radio luminosity has a substantial dispersion. Fu rthermore, it has l ong been suggested 
that the synchrotron rad iation from low SFR galaxies is som ewhat suppres sed (|Klein et al.l 11984 iChi fc Wolfendalel 
119901 : [Price &; D uric 1992), although this view is controversial dCond on 1992). A compar ison of the local SFR fu nction 
with the local radio luminosity function at 1.4 GHz (dominated by synchrotron emission: iMauch ~&T Sadler 2007) shows 
that the latter is flatter than the former at low luminosities. If the corresponding portion of the SFR function was 
inferred from far-IR data one could argue that it is overestimated since a fraction of the far-IR emission could be due 
to dust heated by evolved stars. But in fact it was derived from Ha measurements and therefore it is not liable to 
such a problem. To reconcile the SF R. functio n with the radio luminosity function we need to assume a non-linear 
A S ync-SFR relation. Following iMassardi et al.l (2 0101 ) we adopt a relation of the form: 


to 


where L sync is given by eq. © and f3 = 3. This equation converges to eq. ([3]) for L sync L* jSync = 0.886L sync (SFR = 
1 M 0 yr -1 ). At 1.4GHz, T*, syn c — 10 28 ergs -1 Hz -1 . Applying this relation to the redshift dependent SFR functions 
presented in Sect. [2] (see Fig. QJ and further adopting for the L ra di 0 -SFR relation a dispersion of cr ra dio = 0.4 dex we 


( A* 5 sync / Async) 


= to = 


^*,sync 


(^★,sync/-^sync) “t - 
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get agreement with the iMauch fc Sadleil (12007!') luminosity function (see Fig. 0 as well as with the sub-mJy counts at 
1.4 GHz, dominated by emission from star-forming galaxies. The model counts, obtained adding those of radio-loud 
AGNs yielded by the iMassardi et ahl (12010H model 7 , are shown in Fig. [3] 


3.2. Calibration of the relation between SFR and free-free emission 
We have rewritten the relation between SFR and free-free emission derived bv lMurphv et al.l (j2012f ) as: 

L„ = 3.75 x 10“ (-|A_) (yTjj) gF.T) exp (-g) 


(5) 


w here T is t he te mperature of the emitting plasma and g(z/, T) is the Gaunt factor for which we adopt the approximation 
of iDrainel rf20rT[) . which improves over earlier approximations 


g(i/,T) = In 



5.960 



7r 




+ exp(l) > , 


( 6 ) 


Zi being the charge of ions and T their temperature. This equation reproduces the iMurohv et al.1 (I2012f) calibration 
at the calibration frequency (33 GHz) for a pure hydrogen plasma (Z. t = 1) and T = 10 4 K. At the other frequencies 
relevant here the differences due to the different approximations for the Gaunt factor are small. We will adopt the 
above values of Z iy and T throughout the paper. 

We have used eq. (0 to convert the redshift-dependent SFR functions to free-free luminosity functions. To compute 
the cou nts of star-f or ming g alaxies at 95 GHz, where an observational estimate has been obtained thanks to the SPT 
survey dMocanu et al.l [20130 . we need to take into account that the measured flux densities include, in addition to 
the free-free emission, contributions from thermal dust and from the synchrotron emission (see Fig. 0. All these 
contributions were taken into account using the Tff-SFR and L sync -SFR relations given above plus the Lir-SFR 
relation appropriate for each galaxy population: eq. 0 for late-type galaxies and 


L m = 3.1xl0“(-?f2_) 

VMfh vr i / 


Mq yr- 


erg s 


(7) 


for pro to-spheroid al galaxies dCai et al.ll2013 ). This calibration is essentially identical to the empirical relation derived 
bv IMurohv et al.l ( 2012 b To get down to any chosen frequency we used the SEDs for the three populations of star¬ 
forming galaxies (“cold” and “warm” late-type galaxies, and high-z proto-spheroidal galaxies) considered bv lCai et al.1 
(2013|), adding the free-free and synchrotron emissions given by eqs. © and 0. 

The model counts turned out to be substantially below those estimated bv iMocanu et al.l (1201311 . The statistical 
approach adopted by these authors, however, did not appear completely reliable. They used the a^o spectral indices to 
classify their sources, taking ct^o = 1-5 as the boundary between dust-dominated galaxies ( 0:220 > 1-5) and radio-loud 
AGNs (synchrotron-dominated, in their terminology). They computed the probability that their posterior value for 
0220 was greater than the boundary value, and interpreted it as the probability for a source to be dust-dominated. The 
95 GHz differential counts were calculated as the sum of the probabilities P(a 22 o > 1-5) that the sources in a given 
flux density (SgsGHz) range are dusty. Since at 95 GHz the dusty galaxies are a small fraction of detected sources, this 
approach may easily overestimate the counts because, due to the large uncertainties on ct^o, many radio loud AGNs 
have a non-zero probability of being dust-dominated. The sum of these low probabilities may give a count comparable 
to, or even larger than the count of the rare truly dusty galaxies. 

This argument has motivated us to re-estimate the 95 GHz counts of dusty galaxies using a safer, although lengthier, 
approach. There are 406 95 GHz sources brighter than 12.6 mJy, the 95% completeness limit of the survey. We have 
checked the SEDs of all of them, collecting the photometric data available in the NASA/IPAC Extragalactic Database 
(NED) plus the Wide-held Infr ared Survey Explorer (WISE: iWright et al.l 1201(1) A11WISE 8 catalogue, the Planck 
Catalogue of Compact Sources (|Planck C ollabor ation XXVTIll [20 1 ill and the IRAS catalogue. We found that only 4 
sources, all with P(a > 1.5) = 1, are indeed dusty galaxies. They are listed in Tableland their SEDs are shown 
in Fig. 0 A fifth source with P(a 220 > 1.5) ~ 1, SJ043643 — 6204.4, is only 0.368 arcmin away from the very bright 
red giant Mira variable star R Doradus. Its SPT photometry may be strongly contaminated by emission from dust 
shrouds that often form arou nd this kind_of stars. For this reason we have discarded this source. 

For comparison, the lMocanu et al.l (|2013H approach yields ~ 12 dusty galaxies with S 95 GHZ > 12.6 mJy. A substantial 
contribution to the difference with our result comes from sources with P(a 220 > 1.5) <C 1. The re-assessed integral 
counts are shown in Fig.0 It is worth mentioning that 3 out of the 4 confirmed dusty galaxies have low frequency radio 
luminosity in excess of that expected from the mean relation between IR and 1.4 GHz luminosities. The difference is 
about twice the dispersion of the relation between th e syn ch rotron luminosity and the SFR, yet within the range of 
values for star-forming galaxies according to lYun. Reddv fc Condonl (1200 lh . It may thus be not statistically significant, 


7 Tabulations of the radio-loud AGN counts yielded by the model at several frequencies are publicly available at 
http://wl.ira.inaf.it/rstools/srccnt/srccnt_tables.html . 

0 http://wise2.ipac.caltech.edu/docs/release/allwise/ 
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especially taking into account a possible selection effect: the higher synchrotron luminosity helps bringing the sources 
above the detection limit at 95 GHz. The parameter q defined by 


q = log 


' FIR/3.75 x 10 12 Hz ' 
_Si.4GHz(Wm -2 Hz -1 ) 


( 8 ) 


where 

FIR(Wm- 2 ) = 1.26 x 1(T 14 [2.58S 60Aim (Jy) + Sioo„m(Jy)], (9) 

has values ranging from 1.85 to 2.06 for the aforementioned galax ies, som ewhat smaller than the mean value of q 
for star-forming galaxies (q = 2.34 with a dispersion a = 0.26: lYun. Reddv fc Condor] 1200 lib yet larger than the 
value gii m = 1.64 proposed bv lYun. Reddv fc Condonl (|2001h as the boundary between starburst-dominated and AGN- 
dominated radio emission. In any case, we cannot rule out alternative explanations such as an additional contribution 
from a weak active nucleus or excess radio emission arising from radio continuum brid ges and tidal tails not assoc iated 
with star formation, similar to what is observed for so-called “taffy” galaxies dCondon et al.ll2002t lMurphvll201§l . 

The SEDs of all the other sources are consistent with the 95 GHz flux density being dominated by synchrotron 
emission and inconsistent with them being dusty galaxies. As illustrated by Fig. [3 this is the case also for the other 4 
sources with P(a 22 o > 1-5) > 0.6, including the one with Pietro > 1.5) = 0.91. The radio-source nature of all objects 
with lower values of P(a 220 > 1.5) is even clearer. 

Figures [6] and [8] show that we o btain good agr e emen t with the new determination of SPT counts of dusty galaxies 
at 95 GHz as well as with those bv lMocanu et al.l ( 12013 ) at 150 and 220 GHz. While the 95 GHz counts are dominated 
by nearby galaxies, those at the higher frequencies include a contribution from high-z strongly lensed proto-spheroidal 
galaxies; their contribution increases with decreasing flux density (in the ranges covered by the SPT surveys) to the 
point of becoming dominant at the SPT detection limits. 

The model has been further tested at 4.8, 8.4, 15 and 30 GHz (Fig. 0. Not surprisingly, the free-free contribu¬ 
tion becomes increasingly important at higher and higher frequencies. While at 1.4 GHz the dominant emission is 
synchrotron, at 8.4 GHz the free-free emission takes over at sub-mJy levels. 


4. PREDICTIONS FOR SURVEYS WITH THE SQUARE KILOMETER ARRAY (SKA) AND ITS PRECURSORS 

4.1. Number counts and redshift distributions 

The plans for SKA Phase 1 (SKA1-MID) will likely include three continuum surveys (in band 2), all at ~ 1-1.4 
GHz: one over 1000 deg 2 with rms ~ 1/iJy/beam, a second one over 30 deg 2 with rms ~ 0.2 /iJy/beam, and a third 
over 1 deg 2 with rms ~ 50nJy/beam dPrandoni fe Sevmourll2014D . In addition there will possibly be a deep survey in 
band 5, at ~ 10 GHz, with rms ~ 0.3 /iJy/beam over ldeg 2 , and an ultra-deep survey again in band 5, at ~ 10 GHz 
with rms ~ 0.03 /zJy/beam over ~ 30arcmin 2 (see also lMurnhv et ahll201 4). 

Deep surveys are also pla nned w ith the SKA precursors, like the Australian SKA Pathfinder (ASKAP) and the 
South African MeerKAT ilNorris et al.ll2013 ). The Evolutionary Map of the Universe (EMU) is a project that will use 
ASKAP to make a deep (5er limit of 50 /iJy beam -1 ) radio continuum survey of the entire southern sky, extending 
as far north as +30°, i.e. covering about 30,000deg 2 at ~1.4GHz. The MeerKAT International GigaHertz Tiered 
Extragalactic Exploration (MIGHTEE) Tier 2 will exploit the MeerKAT Phase 2 to observe over 35 deg 2 down to a 
5cr limit of 5/zJybeam - , again at ~ 1.4 GHz. 

The SKAl-MID, EMU and MIGHTEE 5 a flux density limits for unresolved sources are indicated by vertical lines 
in Fig. [3] The contributions of the various dusty galaxy populations are shown by the different lines. The “bump” 
at tens of ^rJy’s is due to late-type galaxies at z ~ 1-1.5. Proto-spheroidal galaxies at higher z become increasingly 
important at lower flux densities, down to a few hundred nJy’s. 

The predicted redshift distributions for surveys at the flux density limits mentioned above are shown in Figs. 1101 
and QH The surface density of galaxies at z > 6 increases rapidly with decreasing flux density. For the 5 a detection 
limits of the 1.4GHz surveys mentioned above it is of ~ 10 -2 , 6.4, 148 and 1013 deg~ 2 for Si. 4 GHz > 50, 5, 1 and 
0.25 ^Jy respectively, implying that the 3 SKAl-MID surveys are expected to detect ~ 6,400, 4,440 and 1013 z > 6 
galaxies, respectively, while the MIGHTEE and the EMU surveys should detect about 220 and 300 galaxies at these 
redshifts, respectively. Figure [TO] shows that the SKAl-MID surveys may reach even higher redshifts, allowing us to 
get a glimpse of the cosmic SFR across the re-ionization epoch, overcoming limitations by dust extinction. We expect 
the detection of about 120, 163 and 58 galaxies at z > 8 by the 3 SKAl-MID surveys down to Si. 4 GHz = 5, 1 and 
0.25 /iJy respectively. 

As illustrated by Fig. QTj the planned band 5 surveys will also reach very high redshifts (see also iMurohv et all 
120141 ). We predict the detection of ~ 16.5 and ~ 5.5 z > 6 galaxies with the deep and the ultra-deep survey, 
respectively. These estimates include the contribution of dust emission that is significant for the hig hest re dshift 
galaxies ( see Fig. m Howeve r, dust emission is likely under-abundant in most galaxies at z > 6 (jMunoz fc Loebl2Qlll : 
iCai et aPl20l4TCapak et al.ll2015D because there is probably not enough time for a substantial dust enrichment of the 
galaxy interstellar medium. However, g alaxies with high SFRs and substantial dust extinctio n have been reporte d 
(|Riechers et al. 2013; Watson et al. 2015), and the timescale for dust enrichment is poorly known (|Mancini et all2015j) . 
SKAl-MID data will provide key information also in this respect. 
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In Figure [T3] we compare the SKA1-MID potential in measuring the evolution of the cosmic SFR to the results of 
Herschel and of the deepest UV and Ha surveys. The planned SKA1-MID surveys can detect galaxies with SFRs 
ranging from less than one Mq /yr for z < 1 to less than hundred M 0 /yr for z < 7. This is an improvement by more 
than two orders of magnitude compared to Herschel results and will allow us to determine the dust-enshrouded star 
formation history up to much higher redshifts. The dee pest rest frame UV s urveys do better, but only for dust-free 
galaxies. The deepest ALMA maps at 1.1 and 1.3 mm (iCarniani et al-lfeOl.T ) reach lower SFRs than SKA1-MID at 
z > 3, thanks to the negative K-correction. We note that ALMA is not designed to be a pure survey instrument and 
has a very small field of view (the Full Width at Half Maximum of the primary beam is 22" at 300 GHz); however it 
is interesting to compare its capabilities with the SKA (see Figure fl2l). 

4.2. Strongly lensed galaxies 

Figures [IU1 and fill also show that a substantial fraction of the highest redshift galaxies that should be detected by 
the foreseen surveys are strongly lensed (magnification fi > 2) protospheroids. The strongly lensed fraction at high 
z i ncreases with increasing flux density limit, while the total counts rapidly decrease. Adopting the SISSA profile 
(lLapi et al.l i2012h for the lens galaxies and an upper limit /i max = 30 to the gravitational amplification, the model 
yields 1195, 432, 101 and 7.6 strongly lensed galaxies per square degree brighter than the SKAl-MID (band 2), the 
MIGHTEE and the EMU 5 cr detection limits of 0.25, 1, 5 and 50 /iJy respectively (see Fig. fTJl) . They are approximately 
the 0.5-0.6% of the total number of galaxies at these flux density limits (to be precise, th e strong ly lensed fra ctions 
are 0.62%, 0.52%, 0.47% and 0.63%, respectively). The model with jt max = 30 was shown bv lBonato et al.1 (12014 their 
Fig. 2) to provide a good fit to the SPT counts of dusty galaxies without an IRAS counterpart, interpreted as high-z 
strongly lensed galaxies. In the case of band 5 SKAl-MID surveys we predict, over the survey areas, cs 4 and ~ 70 
strongly lensed galaxies brighter than 0.15 and 1.5 /zJy respectively. They are about the 0.36% of the total number of 
detected galaxies at both flux density limits. 

To estimate the number of galaxies for which SKAl-MID can provide direct evidence of strong lensing by detecting 
at least two images we have turned to the classical Singular Isothermal Sphere (SIS) profile for lens galaxies, since 
the application of the SISSA model for this purpose is much more cumbersome. Details on the calculation are given 
in the Appendix [B] We find that two images will be detected for about 30% of strongly lensed galaxies brighter than 
0.25/iJy at 1.4 GHz; this fraction decreases to about 20% at a flux limit of 5 /iJy. This means that high-resolution 
radio observations alone will allow us to directly discern strong lensing for 20-30% of the strongly lensed sources. Note 
that the detection of multiple images can happen thanks to the high spatial resolution of SKAl-MID surveys. It will 
not be generally possible with the lower resolution of EMU/ASKAP and MIGHTEE/MeerKAT surveys. 

5. SYNERGIES WITH OPTICAL/NEAR-IR SURVEYS. 

A full scientific exploitation of the data provided by the SKA surveys requires redshift determinations of the detected 
galaxies, especially of the most distant ones which carry information on the earliest stages of galaxy formation. Several 
deep optical/near-IR imaging surveys covering substantial sky areas are ongoing or foreseen. 

In Fig. [14] we compare the minimum SFR of a dust-obscured galaxy (whose SED is shown in the right-hand panel of 
Fig, m detectable, as a function of redshift, by the deepest SKAl-MID survey with those detectable by the following 
surveys: 

• The Euclid deep survey covering an area of 40 deg 2 distributed over two fields close to the North and South 
ecliptic poles to 5er magnitude limits for point source detection of 26 mag in the Y, H and J bands and of 
27.5 mag in the VIS (jLaureijs et al.ll2014 ). 

• The Subaru Hyper Suprime-Cam (HSC) ultra-deep survey 9 dMivazaki et al.l[2012f ) covering two fields of 1.8 deg 2 
each to 5cr point source AB magnitude limits of g = 28.1, r = 27.7, i = 27.4, z = 26.8 and y = 26.3. 

• The 4 Large Synoptic Survey Telescope (LSST: ITvezic et all[2008 1 deep drilling fields of 9.6 deg 2 each, aimed at 
reaching 5 a point source AB magnitude limits of 28 mag in the u, g, r, i and z bands and of 26.8 in the y band 
(see Sect. 9.8 of the LSST Science Book 10 ). 

• The Deep-Wide Survey fDWS: iWindhorst et al.ll2009i) to be carried out with NIRCam on the James Webb Space 
Telescope (JWST). The DWS will image the sky with up to 8 filters. A severe limitation to the exploitation 
of the SKA1-MID/DWS synergy is the smallness of the DWS area, lOOarcmin 2 , i.e. a factor of 1/36 ~ 0.028 
smaller than the area of the deepest SKAl-MID survey (ldeg 2 ). We expect that the SKAl-MID will detect, 
within the DWS area, ~ 6, 200 galaxies at z < 6 but only ~ 14 at z > 6. 

The depth of the 4 LSST deep drilling fields is enough to detect all SKAl-MID dust-obscured galaxies in at least 3 
bands up to z ~ 6 and in at least one band up to z ~ 8. Thus, photometric redshift estimates can be obtained up to 
z ~ 6. With the obvious exception of the JWST/DWS the other surveys are somewhat less expedient. In the case 

9 www.subarutelescope.org/Projects/HSC/surveyplan.html 

10 www.lsst.org/files/docs/sciencebook/SH_whole.pdf 


















C. Mancuso et al. 


of the UV-bright SED (left-hand panel of Fig. fl5l) all galaxies detected by the deepest SKA1-MID survey are easily 
detected by the above optical/near-IR surveys in all bands, so that accurate photometric redshift estimates can be 
obtained. Since, as argued above, strong dust obscuration appears to affect only a minority of faint z > 6 galaxies, 
it may be expected that photometric redshift estimates will be missing only for a small fraction of SKA1-MID high -z 
galaxies. 

One may then wonder what the SKA1-MID adds to the information provided by deepest optical/near-IR surveys. 
Apart from the fact that the SFR inferred from SKA data are immune to the effect of dust obscuration, even in 
the case of no obscuration the SKA measurements are complementary to the optical/near-IR ones. As noted in sub¬ 
sect. 13.11 while the rest frame UV emission measures the instantaneous SFR, the synchrotron emission measures the 
SFR ~ 10 7 yr earlier. Thus, in the absence of dust obscuration, the comparison of the SFR inferred from (rest-frame) 
UV ph otometry with that from low-frequency radio data, dominated by synchrotron, yields a measure of the starburst 
age dClemens et al.lf2008f ). Starbursts younger than ~ 10 7 yr have a synchrotron emission decreasing with the starburst 
age, so that the inferred SFR is lower and lower than that inferred from UV data. By the same token, the comparison 
of SFRs derived from UV and low-frequency radio measurements is informative on variations of the SFR on timescales 
of ~ 10' yr. A note of caution is in order however. This comparison could not work at very high redshifts where the 
syn chrotron emission may be highly suppressed due to th e strong inverse Compton losses off the CMB (see sub-sect. 13.II 
and lNorris et alJl2013HCarilli et al .1120081 : lMurphvll2009l) . so that even at relatively low frequencies the radio emission 
may not be synchrotron dominated. 

A rough redshift estimate for SKA1-MID galaxies missing better information can be obtained from the expected 
correlation of the radio spectral index with redshift. As 2 increases the spectral index between the frequencies v\ and V 2 
in the observer frame, a = \og{S{v\) / S{v 2 )) / log( 1 ^ 2 /^ 1 ), flattens with increasing z, due to the increasing contribution 
of free-free compared to synchrotron. As an example, Fig. Qf)] shows the redshift dependence of the spectral index 
between 1.4 and 4.8 GHz. The accurate photometric redshifts that can be provided by the optical/near-IR surveys 
will allow us to better assess this relation. 


6. CONCLUSIONS 

We have produced detailed predictions for the counts and the redshift distributions for planned surveys with the 
SKA1-MID and its precursors, taking into account the contributions of the different star-forming populations (late- 
type - normal and starburst - and proto-spheroidal galaxies), allowing for the effect of gravitational lensing. As 
mentioned in Sect. [2] the star formation is still ongoing in late late-type galaxies, while it mostly occurred at z > 1.5 
in proto-spheroidal galaxies. Thus the dominant star forming population gradually changes around z ~ 1.5. The 
predictions rely on the redshift evolution of the SFR function inferred from the wealth of data provided by surveys at 
far-IR to mm wavevelengths as well as in the UV and in Ha. These data sets complement each other by probing the 
dust obscured and unobscured star formation, respective ly._ 

In p ractice, we have exploited the models bv lCai et all ( 2013 ?) for the dust obscured star formation and bv lCai et al.l 
( 2014 ) for the total SFR function, including the unobscured one, at z > 2. The missing contribution of the unobscured 
star formation at z < 2 has been added using a phenomenological approach. The model extended in this way has been 
further successfully tested against observational determinations at several redshifts of the Ha luminosity function, 
which is an independent measure of the SFR function. 

The SFR functions were converted into radio luminosity functions takin g into account both the synchrotron and 
the fre e-free emissions, using the calibrations bv lMurnhv et al.1 (|201lLl2012| ). A comparison with the iMauch fc Sadler! 
(|2007f ) local luminosity function for star forming galaxies at 1.4 GHz showed evidence of a deviation from a linear 
relation between the SFR and the synchrotron luminosity, in the sense that low-SFR galaxies are under-luminous. 
This deviation was taken into account by means of a simple analytic formula. Good agreement with the bright portion 
of the luminosity function is achieved adopting a dispersio n of cr ra dio = 0.4 dex for the T r adio _ SFR, relation. 

This procedure, combined with the iMassardi et al. ( 2010) model for radio-loud AGNs, allowed us to reproduce the 
available source counts in the range 1.4-30 GHz, but under-predict ed the published SPT counts of dusty galaxies at 
95 GHz ([Mocanu et al.112013 ). However a careful re-analysis of the iMocanu et all ( 2013 ) sample has shown that the 
surface density of dusty galaxies was appreciably overestimated. The revised counts are in very good agreement with 
model predictions, although the uncertainties are large, due to the poor statistics. Our re-analysis also showed that 
the main contribution to the observed 95 GHz flux density of local galaxies comes from free-free emission, not from 
thermal dust as implied by some models. 

With the radio luminosity function in hand, radio source counts and redshift distributions at the relevant frequencies 
and detection limits were straightforwardly computed. The highest redshift tails of the distributions at the detection 
limits of planned SKA1-MID surveys were found to include a substantial fraction of strongly lensed galaxies. We 
predict that a survey down to 0.25 /rJy will detect about 1195 strongly lensed galaxies per square degree, at redshifts 
of up to 10. For about 30% of them the SKAl-MID can detect 2 images, thus providing direct evidence of strong 
lensing. The integral counts of strongly lensed galaxies are roughly proportional to 5 ,_0 ' 5 in the flux density range 
0.01 /iJy < 5*1. 4GHz < 1 juJy, implying that, at fixed observing time, for the purpose of detecting these objects it is 
more convenient to widen the covered area rather than going deeper. The counts substantially steepen above 10 /iJy. 
Yet the EMU survey covering 30,000 deg 2 to a 50 ^Jy 5 er detection limit is expected to detect ~ 2 x 10 4 strongly lensed 
galaxies. 

The redshift distributions also show that the SKAl-MID surveys, and to some extent surveys with the SKA precur- 
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sors, ASKAP and MeerKAT, will allow us to probe the star formation history through the re-ionization epoch without 
being affected by dust extinction. The deepest SKA1-MID surveys, down to 0.25 /iJy, will extend the determination 
of the SFR function by about two orders of magnitude compared to Herschel surveys up to the highest redshifts. 

Finally, we have discussed the synergies between the deepest SKA1 -MID survey and deep ongoing or planned 
optical/near-IR surveys. We find that the LSST deep drilling fields reach faint enough magnitudes to detect essentially 
all SKAI-MID galaxies in at least 3 bands, allowing reliable photometric redshift estimates. Exceptions may be heavily 
obscured z > 6 galaxies that are however expected to be rare. The JWST/DWS can detect also those in at least 6 
filters, but its small area is a serious limitation. 
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APPENDIX 

EVOLUTION OF THE UV AND Ha LUMINOSITY FUNCTIONS OF LATE-TYPE GALAXIES 


For the UV ( 
EH et, all (2013 1 


A = 1500 A) luminosity functions of late-type galaxies we adopted the functional form chosen by 
for the IR luminosity function of these galaxies: 


$(logTi5oo,z)d logLisoo 


/AisooN 1 “ 

log“(l + Li50o/A*)l 

{ L* ) X6XP 

2 <t 2 


log Ui500- 


(Al) 


We obtained a sufficiently good description of the data using a simple pure luminosity evolution model \L*(z) = 
Lq(1 + z) ql up to 2 = 1, see Fig. [17]. The best-fit values of the parameters are log[4>g (dex -1 Mpc -3 )] = — 2.150±0.095, 
log(Lo/T 0 ) = 9.436 ± 0.119, a = 1.47 7 ± 0.050, a = 0.326 = 1= 0.035, and ar, = 2 .025 ± 0.063. The free par ameters are 
constrained using the data at z ^ 1 bv lWvder et, all (12005H . iTrever et all (1200 5H. a n dIArnouts et all ( 2005 ?). In Fig. [TT] 
the thin solid black lines show, for comparison, the Schechter fits bv fArnouts et, all (120051) . The dot-dashed orange lines 
sho w the contri butions of proto-spheroidal galaxies derived adopting Ai 50 o = Ai 350 -fc(1500 A)/fc(1350 A) = 0.936Ai 350 
dCalzetti et al.ll2000l). 

As argued bv iHopkins et al.l (120011). dust r e ddeni ng of late-type galaxies appears to be dependent on the SFR. 
Adopting their relation and the lCalzetti et al. (l2000l ) extinction law, the extinction of the Ha line is given by: 


A HQ (logSFR) 


6.536 log 


0.797 log(SFR) +3.834 

2^88 


(A2) 


The corresponding UV extinction at 1500 A, based on the lCalzetti et all (| 2000 D law is 

A 1500 = — -WITTd— ) star A na = 3.107 x 0.44 x Ah« = 1.367AH a; (A3) 

/fc(6563A) E(B — U) gas 

where we have taken into account that the color excess appropriate for the stellar continuum, E(B — VUtan is related 
to that for nebular gas emission lines, E(B — V) gas , by E(B — U) s tar = 0.44£ , (H — U) gas dCalzetti et al.ll2000i) . 

The UV luminosity functions of late- type gala xies h ave b een corrected for the SFR dependent dust extinction and 
converted into SFR functions using the iKennicutt & Evans! ( 2012 ) calibration 

l0g (^) = l0g ( C-O + 9 ' 825 " °- 4Al50 °( log SFR )’ (A4) 

via the equation 


d 2 N(log SFR) _ d 2 N(\ogL 1500 ) rflogLi 50 g 
dlogSFRdV dlogLi 5 oodU d log SFR. 

_ d 2 N(log L 150o) L _ q ^ cb4 150 o 
dIogLisoodU d log SFR 


(A5) 


The results are shown in Fig. |l] As a further test of the model we have used the SFR fun ctions obtained co m bining 
UV and IR data to derive the Ha luminosity functions at several redshifts, adopting the IKennicutt &; Evans] (|2012f ) 
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calibration. To compare them with observations we have applied the SFR dependent dust attenuation of eq. (IA2I) in 
the case of late-type galaxies and the Calzetti relation between Ha and UV extinction in the case of proto-spheroids: 


fc(6563A) E(B — U) gas 3.3258 1 

~ A;(1350 A) E{B-V) ste 1350 ~ 11.039 0.44 


0.685Ai35Q. 


(A6) 


As shown by Fig. 1181 the agreement is generally good, with a hint that also Ha surveys may miss galaxies with the 
highest SFRs, particularly at z ~ 1.5. Figure Q] shows that SFR estimates from Ha luminosities exceed those from UV 
luminosities if both are unco rrected for dust attenuation, consistent with the lower attenuation at the Ha compared 
to UV wavelengths [eq. (1A6|) ]. 


COUNTS OF STRONGLY LENSED SOURCES 

In a Singular Isothermal Sphere (SIS) model the strong lensing regime (i.e. the formation of multiple images) 
corresponds to /z to t > 2, where fit ot is the total amplification, sum of the modules of the amplifications of the two 
images. In fact, in a SIS model, two (and only two) images can form in the strong lensing regime. If we indicate with 
fi + and //_ the values of the amplifications of the two images we have that 


M+ = 


R- =' 


0 + $E 

T~ 

0-O E 


0 


> 0 


< 0 


(Bl) 

(B2) 


where 0 is the angular separation of the background source from the optical axis of the lens system. In the strong 
lensing regime 0 < 6e, 0e being the Einstein radius 


6e = 47T 


m 


Dls 
D S ' 


(B3) 


Here a v is the ID velocity dispersion of the lens, D$ is the angular diameter distance to the source, and Des is the 
angular diameter distance from the lens to the source. From the above equations it follows that 


I | 20 e 

M tot = fi+ + \fi-\ = — 


(B4) 


It is worth noticing that the strong lensing regime implies 0 < d E and, therefore, fi+> 2 and |/i_| > 0. 

The cross-section for strong lensing with total amplification higher than a given value //tot > 2 is 

ao2 

s(> Mtot) = 7T-2^- (B5) 

Mtot 

Similarly, the cross section for strong lensing with the second image (the fainter one) having the module of the 
amplification higher than a given value \fi-\ is 




£(> |li-|) =T——|-pr 

(1 + |M-I) 


(B6) 


Once the lensing cross section is known, the corresponding probability is given bv (in th e “non -overlapping” cross- 
sections approximation, i.e. neglecting lensing by more than one foreground mass; IPeacoc sum 


rM? 


P{> //|z s ) = 


E(> fi\M L , ze, z s ) 


d 2 N h 
dMhdzL 


dM L d^L, 


(B7) 


where 2s and Zl are the redshifts of the background source and of the lens, respectively, while d^NE/dM^dzE is the 
number density of lenses per unit mass and redshift interval. The usual choice for th e mass distrib ution of the lenses 
is the mass function of the dark matter halos derived from N-body simulation ('e.g.. ISheth fc Tormenl[l999i ). 

The halo mass is related to 6e (and therefore to the the lensing cross-section) via the ID velocity dispersion of the 


lens, 


by 


a, 2 = 


GM l 

2I?L,vir 


(B8) 


where Re vir is the virial radius of the lens that we compute assuming a virialization redshift ze vir = 2.5 Isee lLapi et all 

!M3). 

Given the true number density of un-lensed sources per unit logarithmic interval in flux density and per unit interval 
in redshift, (d 2 V/dlogS , d2)T, the corresponding number density of lensed sources is computed as 


rd 2 v(s,2 S )l 

/‘Mm ax 

"d 2 N(S//i,z s y 

dlogS d^s 

L ** Mm in 

dlogS” dzs 


p(fi\z s )dfi, 


(B9) 
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where S is the measured flux density of the sources and p(p\zs)dp is the probability that a source at redshift zs ha s its 
flux density boosted by a factor p, within dp. because of a lensing event. This probability is obtained from eq. (lB7l) 
as p(p\zs) = —dP(> p\zs)/dp. The upper integration limit, p maxi is set by the size of the background source, with 
/-Hot,max = oo for a point source. 

The solid and the da shed red curv es in Fig. [13] have been derived by plugging into eq. (IB7I) the lensing cross- 
section given by eq. (1B5I) and Eq. (IB6I) . respectively, and assuming in both cases /ttot,max = 30 (translating to \p-\ < 


/Rot, max/2 1 — 14). 
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Figure 1. SFR functions at several redshifts. The blue lines show the SFR functions derived from the UV luminosity functions given 
by the phenomenological model for late-type galaxies with and without the extinction correction (dashed and thin solid lines, respec¬ 
tively). The dot-dashed orange lines show the contributions given by the physical model for proto-spheroids with a minimum halo mass 
* Mq . The thick black lines show the sum of the contributions fro m late-type galaxies, after the extinction correction, and of 


2010 

Laoi et al. 20113; [GruDoioni et al. 2013 

Maenelli et al. 2011. 2013]). The blue data points are from UV measurements (Wvder et al.l 

2005 

Arnouts et al. 2005|; Heinis et al. 2013|; 

Sawicki & Thompson 200G; Reddy Steidel 2009; Alavi et al. 2014). The orange data points 


2003yGeac let^LJ 20081: IVillar et al.H2008l : I Haves et al .1120101 : ILv et aO20071 . 1201U : IShim et al.H2009l : ISobral et al .1120091 : IColbert^etaLjl^Ol^: 

Sobral et a . 2013: An ct al. 2014). The conversions of attenuated UV and Ha luminosities to SFRs are done using the lKennicutt & Evans 

(120121) calibrations, i.e., log(SFR/[MQ yr -1 ]) = log(LFuv/R©]) — 9.825 and log(SFR/[M 0 yr -1 ]) = log(LH Q ;/[ er g s_1 ] > ) —41.270, respec¬ 
tively. The conversion of the IR luminosity to SFR was done using eq. CD which takes into account that a substantial contribution to 
dust heating in galaxies with low SFRs comes from old stellar populations. Comparing the SFRs estimated from the UV or from the Ha 
luminosities with those inferred from the IR lum inosities it is appar ent t hat the former are incr easingly underestimated with increasing 
SFR and redshift, consistent with the findings by I Reddy et alj (120121) and IBurgarella et al.l (120131) . 
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Figure 2. Synoptic view of the evolution of the cosmic SFR function as given by the model described in Section \2\ The luminosity 
function at each redshift includes the contributions of normal late-type, starburst and proto-spheroidal galaxies. 
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Figure 3. Euclidean normalized differential counts at 1.4 GHz compared with model predictions. Above ~ 1 mjy the counts are 
dominated by canonical radio-loud active galactic nuclei ( AGNs) w hile star-forming galaxies take over at fainter flux densities. 
The dashed brown line shows the model for radio AGNs bv IMassardi et all (I2010T) . The other dashed lines show the contributions 
of the star-forming populations considered in this paper (unlensed and strongly lensed proto-spheroids and late-type galaxies). 
The two triple-dot dashed li nes show the total synchrotron and free-free emissions from these populations. References for the 
data points can be found in.de. jjOtfi_et_al, 1 201011. Th e cy an and the black polygon s show the ranges of 1.4 GHz counts consistent 
with the P(D) distributions bv lGondon et a .1 (1201211 and IV crnstro m et ahl (|2014f l. respectively. Note that thes e P(D ) analyses 
rule out the sub-mjy counts as high as the IQwen fe Morrisonl (1200811 estimates. As argued by ICondon et al.l (1201211 the high 
counts may result from the large and difficult corrections for the effects of partial resolution; the quoted P(D) data do not need 
significant resolution corrections. The vertical lines correspond, from left to right, to the 5 ct detection limits of the planned 
SKA1-MID surveys to unresolved sources: 0.25/rJy (blue), 1 /rjy (yellow), 5 fj,3y (pink). Given the resolution of these surveys 
(~ 0.5": iPrandoni fc SevmouHl2014 'l it is likely that a fraction of sources (especially bright, nearby galaxies) will be resolved out 
so that the surveys will not be 100% complete at the quoted limits. However, the array configuration proposed, which delivers 
an essentially flat point source sensitivity over a large range of synthesized beams, will overcome this problem. The 5 fi Jy limit 
applies also to the Tier 2 MIGHTEE survey (MeerKAT), while the dark green line shows the EMU (ASKAP) survey 5cr limit. 
Again the limits refer to unresolved sources. The MIGHTEE and EMU surveys have a lower resolution than the SKA1-MID 
ones: 8.3"-3.5" in the case of Tier 2 MIGHTEE (http://www.ast.uct.ac.za/arniston2011/vdheyden.pdf) and 10" in the case of 
EMU (1 Norris et alJlMUl . 
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Figure 4. Synchrotron luminosity function of local star-forming galaxies derived from the local SFR function using the non-linear Z/ ra dio — 
SFR relatio n (cf . eq. [4jl and a dispersion of cr ra( jio — 0.4 dex around the mean relation, compared with the observational determination by 
IMauch & Sadler! (120071) . The luminosity function obtained assuming a simple proportionality between L ra( ji 0 and the SFR [eq. (3}] is also 
shown for comparison (dot-dashed line). The dispersion <7 ra di 0 raises the bright portion of the luminosity function while eq. J4J) flattens 
the faint part. 
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Figure 5. SEDs of 95 GHz selected dusty galaxies in the SPT sample by IMocanu et all (120130 : SPT-SJ041736-6246-8 (top left), SPT- 
SJ041959-5456-2 ftop right), SPT -SJ044540-5914-7 (bottom left), SPT-SJ213629-5433-4 (bottom right). The red curve is a fit made by 
eye summing the ICai et al.l (120130 SED of a cold late-type galaxy (dotted blue line), chosen because the 4 galaxies are all normal disc 
galaxies, with the free-free emission given by eq. © for the same SFR (long-dashed magenta line) and with the synchrotron emission 
matching the low frequency data. For three of the galaxies, such synchrotron emission (short-dashed orange line) is higher than that given 
by eq. 7 represented by the long-dashed green line. The di ffere nce is about twice the dispersion of the relation between the synchrotron 
luminosity and the SFR (sha ded green band; see sub-sect. l3.lL yet within the range of values for star-forming galaxies according to 
lYun, Reddy fc Condonl (120011) . It may thus be not statistically significant, especially taking into account a possible selection effect: the 
higher synchrotron luminosity helps bring the sources above the detection limit at 95 GHz. However, we cannot rule out alternative 
explanations such as an additional contribution from a weak active nucleus or excess radio emission arising f rom radio continuum bridges 
and t idal tails not associated with star formation, similar to what is observed for so-called “taffy” galaxies dCondon et al.lf2002l : IMurphvl 
120131 ). Data points: low frequency radio flux densities from the NED (open circles); SPT flux densities (triangles); Planck flux densities 
(open diamonds); IRAS flux densities (open stars). 
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Figure 6. Our re-assessment of the integral counts of dusty galaxies at v =95 GHz (asterisks with error bars) compared with the expec¬ 
tations from the lCai et al.l 1120131) model, using for each source populations the complete SEDs that include dust, synchrotron and free-free 
emissions. The various lines correspond to the contributions from the different source populations and from the different emission processes, 
identified in the legend inside the figure. 
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Figure 7. SPT sources selected at 95 GHz that have 1 > P(a 220 > 1.5) > 0.6: SPT-SJ043651-5841-1 (top left), SPT-SJ203730-6513-3 (top 
right), SPT-SJ204101-5451-4 (bottom left), SPT-SJ213406-5334-3 (bottom right). Their SEDs are incompatible with those of star-forming 
galaxies, illustrated by the solid red line, and consistent with being radio sources. Data points: low frequency radio flux densities from the 
NED (open circles); SPT flux densities (triangles); WISE flux densities (open squares); 2MASS flux densities (spirals); in violet and dark 
green are IRAS and Planck upper li mits, respective ly. We hav e adopted as Planck upper limits the 90% completeness limits of the Planck 
Catalogue of Compact Sources llPlanck Collaboration XXVIHH201 4?) in the “extragalactic zone”. 
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Figure 8. Euclidean normalized differential counts of dusty g alaxies at 150 a nd at 220 GHz (left and right panel, respectively). The 
asterisks with error bars are the observational determinations bv IMocanu et ahl (120130 . The various lines correspond to the contributions 
from the different emission processes and from the different galaxy populations, identified in the legend inside the panels. The lines 
corresponding to the synchrotron emission are below the lower limit of the y-axis. Note that in the considered flux density ranges the 
counts of star-forming galaxies seen via their radio emission are made by nearby sources and have thus an Euclidean slope. On the contrary, 
counts due to unlensed and lensed proto-spheroidal galaxies, seen via their dust emission, are due to strongly evolving high- 2 : objects. 
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Figure 9. Predicted versus observed counts at 4.8, 8.4 and 30 GHz. Dusty galaxies come up at sub-mjy flux density levels and their 
counts are accounted for by the model. At highe r flux densities the counts are dominated by canon ical, AGN powered radio sources. 
The models shown are from [Massardi et al.l (|20101 1 at 1.4, and 4.8 GHz and from fde Zotti et al.l (|2005ll at 8.4 and 30 GHz. The various 
lines correspond to the contributions from the different source populations (dashed yellow lines: late-type galaxies; dashed blue lines: 
unlensed proto-spheroidal galaxies; dashed green lines: strongly lensed proto-spheroidal galaxies; dashed brown lines: radio-loud AGNs; 
solid magenta line: total) and from the different emission processes (triple-dot dashed lines, red for synchrotron, violet for free-free). The 
solid vertical lines indicate the predicted flux density limits, for unresolved sources, of possible deep (yellow) and ultra-deep (blue) band 5 
surveys with SKA1-MID. 
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Figure 10. Predicted redshift distributions at the 5cr detection limits of the deepest SKA1-MID, of the MIGHTEE (MeerKAT) and of 
the EMU (ASKAP) 1.4 GHz surveys (see text). The lines show the contributions of the various galaxy populations: late-type galaxies 
(yellow); unlensed proto-spheroidal galaxies (blue); strongly lensed proto-spheroidal galaxies (green); the magenta line shows the total. 
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Figure 11. Predicted redshift distributions at the 5(7 detection limits of SKAl-MID 10 GHz surveys (see text). The lines show the 
contributions of the various galaxy populations: late-type galaxies (yellow); unlensed proto-spheroidal galaxies (blue); strongly lensed 
proto-spheroidal galaxies (green); the magenta line shows the total. 
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Figure 12. Minimum SFR detectable by SKA1-MID surveys 1.4 GHz down to 0.25 /xJy (dashed magenta line), 1 ^ujy (dashed bottle green 
line) and 5/iJy (dashed blue line), as a function of z, compared with the minimum SFR detected by Herschel /SPIRE surveys, basically 
down to the confusion limits (dashed-dot-dot-dot lines, green for 250 /mi, orange for 350 fim and red for 500 (i m) and by UV/Ha surveys 
(horizontal cyan line). The latter line corresponds to the average minimum absolute magnitude reached by the deepest surveys, that turns 
out to be almost independent of redshift. Note that, while very deep in terms of SFR, the UV/Ha surveys miss the dust-enshrouded star 
formation. Also shown, for comparison , is t he minimum SFR detectable in the deepest ALMA maps at 1.1 and 1.3 mm available so far 
(dot-dashed brown line; ICarniani et al.H2015l) . 
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Figure 13. Predicted integral number counts of proto-spheroids at 1.4 GHz. The solid black curve represents the un-lensed counts. The 
red curves represent the counts of strongly lensed galaxies as a function of the total (sum of all images) flux density (solid) and of the flux 
density of the fainter of the two images (dot-dashed). The vertical lines represent the three SKA1-MID 5 cr detection limits for unresolved 
sources (i.e 0.25/rJy, 1 /zjy, 5/xJy). 
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mum SFR detectable by the deepest SKA1-MID survey at 1.4 GHz (down to 0.25/xJy, dashed magenta line), compared 
i SFR detectable in each spectral band by the Euclid deep survey, by the Subaru HS C u ltra-deep survey, by the LSST 
and by the JWST DWS for the obscured SED shown in the right-hand panel of Fig. 1151 
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Figure 15. Left panel: Detection limits of deep optical/near-IR surveys (left) and of the deepest SKAl-MID survey (magenta symbol 
on the right) compared with the spectral energy distribution of an UV-bright young galaxy with a SFR = 2OM0yr —1 at z = 2 (solid), 
4 (dotted) and 6 (dashed). We have considered the following surveys (from top to bottom): Euclid deep (red), Subaru HSC ultra-deep 
(yellow), LSST deep drilling fields (violet) and JWST (green). The cyan point corre sponds to the flux density limit of the deepest ALMA 
maps at 1.1 and 1.3 m m available so far (dot-dashed brown line: ICarniani et ah 2015). The SED was computed using the GRASIL package 
(Silva et al. 1998) for a Chabricr (2003) Initial Mass Function (IMF), a galaxy age of 10 Myr and a metallicity of 0.1 Zq. Right panel: same 
as the left panel but for an obscured SED with the same IMF and SFR, a galactic age of 1 Gyr and a metallicity of 3 Zq. 
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Figure 16. Spectral index between 1.4 and 4.8 GHz, a ( S u ex i /~“), in the observer frame, resulting from the combination of synchrotron 
and free-free emissions, as a function of 2 . The measured spectral index is thus a rough redshift indicator. Note, however, that the plotted 
curve does not take into account the possible increase of the synchrotron suppression in galaxies with increasing redshift due to inverse 
Compton scattering off the CMB photons (see Sect. This effect could hasten the convergence of the mean ol(z) towards the free-free 
spectral index (a ~ —0.1) and introduce significant scatter especially at high 2 : since this effect is highly dependent on how compact the 
star formation activity is within a system. 
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Figure 17. Comparison of the observed 1500 A luminosity functions at several redshifts with the model described in the text. The 
dashed blue lines correspond to the pure luminosity evolution model for late-type galaxies. The contribution of proto-spheroidal galaxies 
(dot-dashed orange lines) begins to show up at z 1.5 and becomes dominant at higher redshifts. This contribution was obtained from the 
1350 A luminosity functions computed bv ICai et al.l (|2O140 assuming a flat intrinsic UV spectrum [e.g., 1/^(1500 A) = 1/^(1350 A)] and a 
ICalzetti et al.l (1200011 extinction curve, yielding A 1500 = ^1350 x fc(1500 A)//c(1350 A) = O. 936 A 1350 . The data at different UV wavelengths 
have been brought to 1500 A in an a nalogous way. I n som e panels we al s o show, for comparison, th e model byE rnouts et al. (2005, thin 
solid black lines) . The data ar e from lArnouts et a l. (2005); Trever ct al. ( 200 5111 ; jWvder e t al. (2005); Sawicki &; Thompson (2006:], ST06); 

I Reddy &; S teidel (2009], RS09); fSteidel et al.ll (1999): Hcinis et al.l (120131) ; I Alaviet alj (120141) ! 
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Figure 18. Comparison of the observed Ha luminosity functions not corrected for dust attenuation with those yielded by the model, 
attenuated as described in the text. The dashed blue lines are for the UV-bright late-type galaxies and the dot-dashed orange lines for 
proto-spheroids. 
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Galaxy 

redshift 

S 95 GHZ (mjy) 

SJ041736-6246.8 

0.00435 

25.56 

SJ041959-5456.2 

0.005017 

24.87 

SJ044540-5914.7 

0.00444 

21.2 

SJ213629-5433.4 

0.002825 

12.76 


Table 1 

Dusty galaxies in the IMocanu et ali lt2013 i 95 GHz sample. 








